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Unconventional electromagnetic mode in neutral Weyl Semimetals 
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We study light propagation in a neutral Weyl semimetal with the Fermi level lying at the Weyl 
nodes in the weak self-interacting regime. The nontrivial topology induces a screening effect in one 
of the two transverse gauge helds, for which we hnd two branches of attenuated collective excitations. 
In addition to the known topologically gaped photon mode, a novel massless and slightly damped 
excitation appears. Strikingly, at low energies this new excitation has a linear dispersion and it 
propagates with the same velocity than the electrons, while at energies well above the electron-hole 
continuum threshold it behaves as a massive attenuated photon with velocity similar to the speed of 
light in the material. There is a crossover at certain momentum in the direction perpendicular to the 
separation of the Weyl nodes above which the novel gapless mode enters into an overdamped regime. 
Regarding the unscreened gauge field we show that it is also attenuated, which is a non-topological 
property shared by Dirac semimetals as well. 
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I. INTRODUCTION 

Weyl semimetals (WSM) are three dimensional topo¬ 
logical electronic systems whose low energy band struc¬ 
ture is described by pairs of Weyl fermions separated in 
momentum space. They were first proposed to occur in 
a class of iridate material^, but it was only recently that 
they were experimentally realizecPM^. The topological 
properties of WSMs are encoded in a Chern-Simons (CS) 
term that stems from the axial anomaly. The anomaly 
is responsible of a plethora of phenomena like the axial 
magnetic effectPMH^ chiral separation effectl^, chiral mag¬ 
netic wavejm and, more recently, a new mechanism for 
the phonon Hall viscosit}^^. 

The experimental evidence of the WSM phase is based 
on ARPES and transport experiments^, and on optical 
conductivity experimentS^ (optical me asurem ents have 
been done in Dirac semi meta ls as welP^^^. The CS 
term induces birefringenc e!^ ^ I ^ and circular dichroisn^^. 
It has been also proposed that chiral electromagnetic 
waves exist and propagate along domain walls in mag¬ 
netic WSMJI^. The interaction of light with magnons in 
WSMs has been addressed iiM while plasmons in Weyl 
metals have also been studiecP^^^. More recently, a new 
type of Helicon mode under applied magnetic fields has 
been proposed^ at finite doping. Here we address the 
problem of light propagation in a WSM at the neutral¬ 
ity point (i.e., the Fermi level lies at the Weyl nodes) 
by computing the poles of the full photon propagator. A 
similar approach was done for Quantum Electrodynamics 
(QED) under an external magnetic fielcP^. Eor the case 
of the solid state-realized WSM there are two key aspects 
that are central: 1 ) nontrivial topology is present even in 
the absence of an external magnetic field and 2 ) contrary 
to QED there are two velocity scales in WSMs: the Eermi 
velocity of the electrons and the velocity of light in the 
material. We will see how these two features combine so 
that novel collective excitation behavior is realized. We 
will resort ourselves on a RPA treatment of the problem. 
This implies that we will effectively consider a non-self in¬ 


teracting Weyl fermions interacting with an external elec¬ 
tromagnetic field. The validity of a RPA treatment has 
been extensively studied in the literature for two dimen¬ 
sional Dirac materialJ^ and more recently in the g eneral 
case of interacting three dimensional Dirac system^^^^. 
Eollowing these references, some of the results presented 
here might be extended to the self-interacting case under 
particular conditions {1/N expansions and so). 

The rest of the paper is structured as follows. In sec¬ 
tion H we describe the model for WSM that we use and 
set notations. In section HI, we give the neccessary de¬ 
tails for the computation of the polarization function and 
explain it. In section IV, we compute the spectrum of 
the electromagnetic collective modes through the spec¬ 
tral function, and describe the most salient features ob¬ 
served. In section V, we discuss the results and mention 
some future lines to follow from this work. 


II. THE MODEL 

Our starting point is the low energy action of a WSM 
coupled to an electromagnetic field: 

S = y d^a;^'(z7°(ao+ieAo)-w7*(5j + ^^j)+6j7*7®)^', 

( 1 ) 

where v is the Eermi velocity and c is the speed of light 
in the corresponding material. This action describes two 
Weyl fermions separated in momentum space by a vector 
2b/v so that time reversal symmetry is broken. Without 
loss of generality will fix 6 = 6 x 3 . Quantum corrections 
to 0 are well known to generate a nonzero odd part of 
the polarization function of the photon in the form of 
a CS ternP^M^. The coefficient accompanying the CS 
term is ill defined, as it depends on the routing of in¬ 
ternal momenta in the loop integrals as well as on the 
regularization method. In a solid state realization of a 
WSM this coefficient can be fixed by means of the bulk¬ 
boundary correspondence relating the theory in the bulk 
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with the Fermi arcs at the edge^^^^. Another way is 
by considering the 3D band structure of a WSM as a set 
of independent 2D band st ructure s, each of them with a 
well defined Chern numbeJ^ ^ l ^^ l ^^ i 

In addition to the action for the fermionic sector Q, 
the action for the electromagnetic field is: 

S = /d‘xl(.|E|»-!^), (2) 

where e = eoe^ and /i = /io/^r are the dielectric permit¬ 
tivity and magnetic permeability of the material respec¬ 
tively. The electric and magnetic fields are related to the 
gauge field as: 

Ei = - diAo, Bi = Kjkd^AK (3) 

III. POLARIZATION FUNCTION 

In what follows we will obtain the 1-loop order polar¬ 
ization function in terms of A^. The computation has 
been performed and deeply examined by several authors, 
both the parity-everP^J® and parity-o d(p8H32] parts of the 
polarization function. We can write it as: 

= + (4) 

where IIq^ is the contribution at zero order in 6, and 11^^ 
is the first order correction. IiJSH 

it was shown that the 
polarization function is analytic in 6, and the correction 
to the even part is second order in 6, breaks gauge invari¬ 
ance and gives a mass to the photon. Imposing gauge 
invariance by using a regularization which preserves it, 
as dimensional or Pa nli-Vi llars regularizations, one ob¬ 
tains a zero correctioiff^ ^. R egarding the odd part, it 
is known to be linear in pHESl and no corrections appear 
beyond the first order. Hence eq. Q is exact, i.e. there 
are not higher order corrections in b. 

The computation of Hq^ is standard. Using dimen¬ 
sional regularization and in the minimal subtraction 
renormalization scheme: 

n(pVMo) = ^(31og(-pVMo) - 5), (6) 

where the four momentum is = (uj^vp), and where 
the space components of Hq^ have a relative weight v/c 
with respect to the time components, as can be read from 
eq. Q. From now on we will redefine p as vp p so 
that ^ — |pp. The parameter /io is an energy scale 

that appears in the process of dimensional regularization 
(see for exampl^^. The appearance of this parameter 
reflects the formal absence of any characteristic scale in 
the linear electronic spectrum. It could be fixed by either 
the experiment or by invoking a standard renormaliza¬ 
tion program (here we are dealing with a modification of 


standard quantum electrodynamics, that is well known 
to be renormalizable) or related to the lattice cutoff of 
the underlying bandstructure by considering a full lattice 
model of the Weyl semimetal, keeping in mind that there 
are several materials with very different lattice structures 
proposed to be Weyl semimetals. For typical values of 
lattice spacings of the order of Angstroms, /io can be es¬ 
timated to be /io ^ 10 eV, however, the reader should 
keep in mind that the low energy/long wavelength as¬ 
pects of the results presented here do not depend on the 
particular value of /io as long as all the frequencies and 
momenta involved in the problem are presented in a di¬ 
mensionless fashion (that is, divided by /io) no matter 
which regularization scheme is used. 

Note that the function n(p^//io) in eq. ® has a 
constant imaginary part for > 0 which defines the 
electron-hole continuum threshold. A normal procedure 
would be to take the local limit oj/\p\ ^ 1 in Hq^. As 
we know for graphene, at the neutrality point and zero 
temperature no well defined plasmon mode exists within 
the random phase approximation (RPA). However, it was 
shown that a plasmon pole appears beyond RPA when 
the momentum dependence in the polarization tensor 
is maintained^, which highlights the importance of not 
considering the local limit in some circumstances. Here 
we will limit ourselves to the RPA approximation but we 
will keep the momentum dependence, and we will see that 
staying away from the local limit is crucial for capturing 
the physics near the electron-hole continuum threshold. 

nr is precisely the CS term. We invok e the boundary- 
bulk correspondence to fix its valu j^^ l ^^ . On the Aq = 0 
gauge nr has only one independent component: 

The presence of this non-diagonal term in the polariza¬ 
tion tensor induces gyrotropy in the WSM, so that left- 
and right-handed polarized light propagate at different 
speedl^. 

Since we are working under the free electron approx¬ 
imation, eq.Q is exact. In an interacting theory, H^^ 
acquires dynamically generated corrections( se^^ for an 
explicit discussion when interactions at constant external 
magnetic field and finite densities are considered, oi^lfor 
corrections in the strong coupling regime). 


IV. COLLECTIVE MODES. 

The inverse full photon propagator can be written as: 

- nr, (8) 

where the inverse free propagator can be obtained from 
eqs. ( |2|3[ ). To obtain the collective excitations we 
look for the poles of the full propagator. The equation 
det[{A^^)~^] = 0 on the Aq = 0 gauge reads: 
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(iX (i-fln(pVMo))(^^-/3r(pVMo)|pP) x 

X ('(i-5n(pVMo))(‘^^-/3r(pVMo)|p|^)=o> (9) 

\ V / \ / flQ flQ /Iq / 


where |p_Lp = Pi + P 2 represents the momentum in the 
plane perpendicular to b, g = /{iirveoer) is the fine 

structure constant, j3 = (?jv^ gives the squared ratio 
between the speed of light and the Fermi velocity, and: 


nv^ipl) 


i-fn(pVM^) 

\-g^{p^lpl)' 


( 10 ) 


J{x,y,z) = 


X 


6 ^ 

TT^ 1 — < 711 ( 0 : — y — z) 

X - j3T{x - y - z)y 
X - l3T{x - y - z){y + z)' 


( 11 ) 


We will take the Fermi velocity of the WSM similar to 
that of graphene {v = lO^m/s) and assume a coupling 
constant g = which is reached with a relative permit¬ 
tivity of the material of ~ 2.3. The speed of light of 
the WSM is given by c = Cr^acuumj (assuming a rel¬ 
ative permeability = 1 ) and the value of the squared 
ratio between the speed of light and the Fermi velocity is 
(3 ~ 3.9 X 10^. Since we are not taking into account inter¬ 
actions between Weyl electrons, a value of the effective 
coupling ^ = 1 is not problematic, and this precise value 
has been chosen for convenience. In realistic solid-state 
WSM, however, it is not expected to have these values 
for vf and/or larger values of e^, resulting effectively in 
an smaller value of g. 

There is a zero mode = 0 corresponding to the un¬ 
physical longitudinal gauge field, arising from the resid¬ 
ual gauge freedom, as x4o = 0 does not completely fix the 
gauge. For the remaining part, there are two physical 
transverse gauge fields. The dispersion relation of the 
unscreened field can be obtained from: 


(i - fln(pVMo)) - /?r(pVMo)|pp) = 0. (12) 

Instead of directly solving the equation, we will define a 
spectral function for this field: 


. = _2^ _ ^ _ . 

[ (1 - 9^{p‘^IpI)) (^2 - (iVip^/pDW) 

(13) 

The spectral function is plotted in Fig. [^a). We see that 
there is a peak with finite width, which corresponds to a 
damped mode, so that the delta function corresponding 
to a free photon is replaced in a WSM by an attenuated 
excitation. Hence our first finding is that light is atten¬ 
uated when propagating through a WSM. This behavior 



FIG. 1: (Color online) (a) Spectral function p for the un¬ 
screened gauge field as a function of the frequency uj and 
momentum |p| in units of po. (b) Spectral function p for the 
screened gauge field as a function of uj and \p 3 \ in units of 
/io, for zero transverse momentum |p_l| = 0. (c,d) Spectral 
function p for the screened gauge field as a function of the 
frequency uj and the transverse momentum | p _ l | in units of 
/io, for (c) IpsI =0 and (d) |p3|/mo = 3 x 10“^. All plots for 
the screened gauge field are presented for iP' / Pq — 3. 


is rooted to the existence of two different velocity scales: 
the value of the Fermi velocity is much smaller than the 
speed of light. This mismatch is encoded in the func¬ 
tion F appearing in eq. ( 12 ), which fulfills lim/ 3^1 F = 1 . 
This means that in the limit where the value of the Fermi 
velocity matches the value of the speed of light in the ma¬ 
terial, we have a non attenuated free photon traveling at 
the speed of light with dispersion relation uj = \p\. Al¬ 
though the gapless nature of WSMs means that there is 
no gap to overcome for electron-hole pair creation, we 
just showed that this is not enough for light attenuation, 
and that damping is generated for a system of gapless 
electrons only when there exist two non-equal competing 
velocity scales. This happens also in Dirac semimetals 
since this property is not related to the CS term. 


The nontrivial topology of the WSM phase shows up 
inducing a screening in the other transverse gauge field. 
The dispersion relation can be obtained from the second 
line of eq. Instead of solving the equation, we will 
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FIG. 2: (Color online) Both plots are presented for a value = 3 and in a logarithmic scale for the horizontal axis, (a) 

Spectral function p for the screened gauge field as a function of the frequency u in units of po for p± = 0 and different values of 
\p 3 1. (b) Spectral function p for the screened gauge field as a function of the frequency uj in units of /xo for ps = 0 and different 
values of \p±\. 


again introduce a spectral function for this field: 


P = 


- 2 ^ 


Mo{(l -5n(pVAio))b^ 


/3r(pVMo)|Pp) - 


Mo ’Mo i 


(14) 


In the analysis following we will consider two separate 
cases: p± = 0 and p± 0. 

i) p± = 0. In Fig. Qb) the spectral function is plotted 
as a function of uj and the momentum jps | in the direction 
of b. Two branches of attenuated collective excitations 
coexist, one of them clearly gaped. The gaped nature of 
the screened field is expected in WSMs, as the CS term 
acts as a topological mass for the photorP^. However we 
see that the second branch corresponds to a new exci¬ 
tation which seems to be gapless and considerably less 
attenuated than its gaped companion. In fact, at a cer¬ 
tain large momentum ps there is a crossover where the 
gaped excitation enters into an overdamped regime, as 
can be seen from the absence of a well defined peak in 
Fig. gb). To highlight the behavior of the novel gap¬ 
less excitation at low frequencies near the electron-hole 
continuum threshold uj = |p|, we plot in Fig. I^a) the 
spectral function as a function of the frequency in a log¬ 
arithmic scale, for different values of ps. We see that 
there is a low energy crossover also for the novel exci¬ 
tation from an underdamped to an overdamped regime 
as the wavelength increases. Overdamping appears for 
frequencies near the electron-hole continuum threshold. 
This suggests a gapless nature of this new excitation with 
dispersion relation uo \p\ at low energies. We see also 
that increasing the momentum ps, the gapless branch de¬ 
viates from the linear dispersion behavior and approaches 
the gaped one. From Fig. Bb) it is clear that the gapless 
branch starts to behave as a massive excitation for large 
enough momentum ps. 

ii) Pj_ ^ 0. Figs. Be ,d) show the spectral function 
as a function of uj and |p^| for the values ps = 0 and 


IpsI/P o = 3 X 10“^ respectively. It can be seen that in¬ 
creasing the perpendicular momentum produces a merg¬ 
ing of the two branches into a single peak. This is a 
result of a second, high energy crossover for the gapless 
excitation which again becomes overdamped, but now at 
big enough p±. For energies near the electron-hole con¬ 
tinuum threshold, we can see in Fig. Bb) that small 
fluctuations in p± drive the gapless excitation across the 
low energy crossover from the overdamped to the un¬ 
derdamped regime, so in a realistic situation when the 
propagation is not perfectly aligned with the b direction 
the gapless branch behaves as a well defined attenuated 
excitation at low energies. 

It is important to note that the behavior near the 
electron-hole continuum would have never been captured 
in the local limit uj/\p\ ^ 1 of the polarization func¬ 
tion. In particular the low energy crossover into the over¬ 
damped regime of the massless branch would have never 
been captured without considering the full momentum 
dependence. 


V. DISCUSSION AND CONCLUSSIONS 

We have studied light propagation in a neutral WSM 
where the Fermi level lies at the Weyl nodes. In this 
situation there are two degrees of freedom corresponding 
to two transverse gauge fields, while the plasmon longi¬ 
tudinal mode that appears at finite chemical potential is 
absent. Besides the generation of a plasmon mode, the 
presence of a plasma frequency would also induce a gap 
to both transverse gauge fields, however at the neutral¬ 
ity point there is only one gaped mode, whose mass is a 
consequence of the nontrivial topology. 

Our first finding is that light is attenuated when prop¬ 
agating through both Weyl and Dirac semimetals. In the 
case of WSMs, nontrivial topology shows up as a screen¬ 
ing effect in one of the two transverse gauge fields, for 
which we find two branches of attenuated collective ex- 
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citations. In addition to the known topologically gaped 
photon mode, a novel massless attenuated excitation ap¬ 
pears, which is our second and most important finding. 
This splitting of a transverse field into two branches was 
obtained before for QED in the presence of an external 
magnetic bleeps, where the ciclotron frequency defines 
the mass threshold for the electron-hole continuum. In 
the case of a WSM there is no mass threshold as nontriv¬ 
ial topology is present in the absence of magnetic fields. 
Furthermore for the case of the solid state-realized WSM 
the Fermi velocity of the electrons is much smaller than 
the speed of light, which translates into a completely 
novel collective behavior of the two branches. 

In particular we find that in a realistic situation of light 
collimated along the direction of b but with small devia¬ 
tions in the perpendicular direction, the massless branch 
shows up as a well defined and damped peak in the spec¬ 
tral function. At low energies its dispersion relation is 
given hy uj v\p\, which defines an excitation traveling 
coherently with the electrons at the Fermi velocity. At 
the other limit of the spectrum {uj ^ v\p\) the gapless 
branch behaves as a massive damped photon with veloc¬ 
ity close to the speed of light. In the situation where light 
propagation is not aligned with b and the momentum in 
the transverse direction is large enough, a crossover ex¬ 
ists and the massless branch enters into an overdamped 
regime, corresponding to a peak with divergent width in 
the spectral function. 

The novel optical excitation could provide a new and 
clean signature for WSMs. It could be detected in exper¬ 
iments as an additional resonance on the transmission 
coefficient of light. There are also indirect proofs, like 


electron energy-loss spectroscopy ffiELS) or subluminal 
Cherenkov radiation spectroscop 3 @^ stemming from the 
small phase velocity of the novel excitation. In addi¬ 
tion, due to the matching of the propagation velocity 
of the novel mode with the Fermi velocity of the elec¬ 
trons, a counter intuitive enhancement of the effective 
dressed electron-photon interaction is expected. This is 
qualitatively different from the effect of the Coulomb in¬ 
teraction enhancement due to the renormalization of the 
Fermi velocity that appears, for instance, in graphen^^. 
It should be noted that, although their effect is negligible 
on the bulk propagation of light, the ID Fermi arc chi¬ 
ral states at the edges of a WSM would give additional 
resonances as surface electromagnetic modes. However, 
the polarization tensor of ID chi ral ele ctrons has already 
been computed in the literatur^^^EH ^j^d therefore the 
Fermi arc contribution can be distinguished from bulk 
propagation. 
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